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Abstract 

The characterization of out-of-local-thermal-equilibrium plasmas requires the use of collisional- 
radiative kinetic equations. This leads to the solution of large linear systems, for which statistical 
treatments such as configuration average may bring considerable simplification. In order to check 
the validity of this procedure, a criterion based on the comparison between a partial-rate systems 
and the Saha-Boltzmann solution is discussed in detail here. Several forms of this criterion are 
discussed. The interest of these variants is that they involve each type of relevant transition 
(collisional or radiative), which allows one to check separately the influence of each of these processes 
on the configuration-average validity. The method is illustrated by a charge-distribution analysis 
in carbon and neon plasmas. Finally, it is demonstrated that when the energy dispersion of every 
populated configuration is smaller than the electron thermal energy, the proposed criterion is 
fulfilled in each of its forms. 

PACS numbers: 52.25-b;52.25.Kn;52.25.Dg 
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I. INTRODUCTION 



Absorption and emission spectra in warm dense plasmas often exhibit broad structures 
described as unresolved transition arrays (UTA) [l|, |2|, |3j, spin-orbit split arrays (SOSA) [4j 
or supertransition arrays (STA) [;5j]. In order to account for huge numbers of transitions, one 
usually resort to statistical treatments based on configuration or superconfiguration average 

ipflnn 

6|, 13, 18|, |9|, H0[ . However, such an average may lack the accuracy required to describe isolated 
levels or lines. A possible improvement lies in the definition of effective temperatures inside 
each subset (configuration or superconfiguration), the detailed level populations inside a 
given subset being derived from the population of this subset and from these temperatures. 



The determination of such temperatures has been discussed in a series of papers 
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13). Alternatively, one may eon.de. hybrid m ode, S involving both fine-strnetnre leve, s and 
configuration-averaged levels [14{ |. 

In addition to the complexity of the spectra, in many situations such as those prevailing 
in laser- or discharge-produced plasmas devised from extreme-UV production 15, Q, \v\ . 
the electron density is not large enough to ensure local thermal equilibrium (LTE) through 
electron-ion collisions. One must then use collisional-radiative models: a set of pertinent 
transition rates is computed and stationary or time-dependant rate equations are then solved. 

Hi; 



OMIC [18|, AT0M3R-0P [19(, AVERROES 
20|, MOST [10], SCAALP [2llL SCROLL Q, Q - - have been developed and 



A series of such models — among which A 
FLYCHK 



benchmarked in the NLTE workshops 22, 



23|. 



While allowing for an accurate description of plasmas in a wide range of conditions, the 
collisional-radiative codes suffer from a practical intractability when large numbers of levels 
- thousands or more — are involved. On must then resort to the above mentioned average 
methods. However, this averaging procedure must be validated. In a preceding paper 



24j . we proposed a criterion based on the solution of a partial-rate system and showed its 



efficiency in checking the validity of the configuration average (CA) procedure. This partial- 
rate system involved electron-impact excitation and ionization, plus the reverse processes. 
Nevertheless, since this criterion only involves collisional rates, it is desirable to implement 
a procedure that includes all types of transition processes. To this respect, a generalization 
of this method is proposed here, where every process is included in the validity check. This 
allows us to diagnose inhomogeneities in all kind of rates. 
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In section [TT1 the general procedure for checking CA validity is presented. Section II II I 
deals with the case where the partial rate system contains radiative processes only. The 
CA-diagnostic method in several of its variants is then illustrated by detailed-versus-CA 
computations in carbon and neon plasmas (section II VI) . Properties of the averaged rate 
system in two particular cases are analyzed in appendices \M and [B] 

II. A VALIDITY CHECK BASED ON DETAILED BALANCE FOR THE 
CONFIGURATION-AVERAGE PROCEDURE 

A. Detailed rate equations 

Let us consider an homogeneous plasma where electrons are assumed to be at thermal 
equilibrium with temperature T e . The ions, with less collisional interactions are out of 
thermodynamic equilibrium when the electron density N e is low. No outer electromagnetic 
field is considered, which amounts to deal with optically thin media or in zero-temperature 
radiation fields. Ion-ion collisions and electron free-free transitions are not included, as they 
weakly affect the population transfers considered here. 

The included processes are radiative deexcitation, collisional ionization, three-body re- 
combination, collisional excitation and deexcitation, radiative recombination, autoionization 
and dielectronic recombination. In the absence of any external field, the radiative deexci- 
tation is unbalanced by the reverse process (photoabsorption) , and radiative recombination 
is not balanced by photoionization. Accordingly, stimulated emission (bound-bound) and 
stimulated radiative recombination (free-bound) are not included either. 

In stationary regimes, plasma properties such as the charge distribution, the average 
internal energy, or the radiative parameters (opacity, emissivity, radiative losses) follow 
from the solution of the so-called collisional-radiative detailed rate equations 




(III) 
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where Rij stands for the sum of all transition rates from i to j. 
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B. Configuration-average procedure 



Assuming the ionic level i (resp. j) belong to configuration a (resp. (3), the average rate 
from a to (3 is defined as 

Rap = — ^2 9i R iV With 9a = ^ 9i (H.2) 

9a 

ISO i£q 

where ^ is the z-level degeneracy, and g a is the configuration degeneracy. If the Ni are 
solutions of the detailed system (111.11) . the total populations N a = J^i^o^i are usually not 
solution of the stationary configuration-averaged rate equations 

</A " -J2 R <*P N <* + J2 R P<* N e = ( IL3 ) 



dt 

P+a P^a 



because the equations fill. 1 j) involve the product RN which is nonlinear. However (see 



Appendix [A} one can derive that the sum of the detailed level populations is indeed a 
solution of the CA equation ( 111.31) in the special case where the transition rates R^ inside a 
pair of given configurations is simply proportional to the final degeneracy gj. 

C. Validity criteria for the configuration average derived from detailed balance 

Intuitively the CA procedure validity should require that the average energy dispersion 
(rms) must be less than the thermal electron energy 

J]iV a AK<T e (II.4) 

a 

where AE a is the rms energy dispersion of the levels belonging to the configuration a, its 
population N a being normalized according to J2 a N a = 1. 

However it has been established in the carbon test case, principally for low T e , that the 



(III. 41) criterion may severely fail 24j. Therefore a different criterion must be elaborated. It 
stems from the detailed-balance (or microreversibility) principle on transitions between any 
pair i, j of (detailed) levels 

AfB^W = N SB R ( g ) (n.5) 
Rf-j being the rate of any transition from i to j and Rjf the rate for the inverse process. In 



equation (1II.5j) . A^ SB is the Saha-Boltzmann (i.e., local thermal equilibrium) population of 



the i level, obeying 



Nf B /N? B = 6^ expi-E^/Te) (II.6) 

9: 



where Eji = Ej — Ei is the transition energy between % (net charge Z ) and j (net charge 

Z*+s), = 2/(N c X^ h ), N e being the electron density and A t h the thermal wavelength 
V(2vrm e T e ) 1 /2. 

In the case where the energy dispersion in every configuration is much smaller than T e , 

AE a < T e for all a (II. 7) 

it can be demonstrated (see Appendix [Bj that microreversibility holds for configurations 
too, i.e., that one has N^ B R^1 = N^ B R^l for every pair of configurations. It implies 
that the configuration population distribution arising from the solution of the averaged 
microreversible rate equations agrees with Saha-Boltzmann law. 



The nature of the processes (p) and (q) involved in equation flll.5j) has not been discussed 



yet. A previous work 24j] had restricted the discussion to the case where Rij are the col- 
lisional rates only. Here we want to compare the microreversibility criteria derived when 
these processes are 



24j . hereafter named 



• collisional excitation and ionization plus inverse processes, as in 
"collisional" case; 

• photoexcitation and photoionization in a Active Planckian field, as discussed in the 
next section ("radiative" case); 

• both above processes plus autoionization and dielectronic recombination ("complete" 
case) . 

Let us notice that this last case is not identical to the usual collisional-radiative case, since 
it involves a fictive Planckian field. 

The interest of these various criteria is that they allow to check the dispersion inside 
a given pair of configurations of the probabilities for each kind of process. This helps in 
determining which process is responsible for a breakdown of CA. For instance, as noted 



before 25], collisional excitation cross sections may be difficult to derive from a simple fit 



formula such as those by Goett et al 26j] , while radiative transition rates are more regular. It 



must be noted that, in the CA case, the computational effort represented by these tests only 



amounts to the solution of one additional linear system for each criterion, since the transition 
rates have already been computed to solve the collisional-radiative problem. When one 
considers CA transitions, these matrix inversions represent little extra computation versus 
the evaluation of the various rates. To sum up, the proposed procedure consists in the 
following steps: computation of the detailed rates for each process with a suitable atomic 
code, followed by the averaging (111.21) ; resolution of the usual CA collisional- radiative system 
(III.3j) ; resolution of one of the modified rate-equation systems as enumerated above and 
comparison with the Saha-Boltzmann solution. This comparison provides a direct indication 
of the dispersion of the transition probabilities within each pair of configurations, and thus 
a "figure of merit" of the CA approximation. 

The above-mentioned modified-rate systems may also be useful in the detailed case as 
check for the matrix-inversion algorithm used: this is illustrated in subsection IIVBI 



III. RATE EQUATIONS INVOLVING RADIATIVE PROCESSES IN A FICTIVE 
PLANCKIAN FIELD 

The purpose of this section is to derive a rate-equation system involving radiative transi- 
tions obeying to the microreversibility principle. This implies that the solution of this system 
must be Saha-Boltzmann in the detailed case (because Saha-Boltzmann solution does satisfy 
the rate equations and because the solution is assumed unique) . In order to allow for detailed 
balance, one must consider here a fictive electromagnetic field with a Planckian distribution 
at a temperature T T = T c . 



A. Definition of the bound-bound radiative rates 



If the levels i and j correspond to the same charge state and Eji > 0, the absorption rate 
from % to j in a fictive field of spectral density u v is related to the spontaneous emission rate 
Aji according to the Einstein relations 

r „, _ 9j A h 2 c 3 _ gj Ajj 

gi 87rEj!j g { e E ^l T - - 1 
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where the density u v is evaluated at the transition frequency Eji/h. In the same way, the 
stimulated emission rate in the u v field is, using again Einstein relations, 

BjiUu = —BijU v = A ^ — -. (III.2) 

J g. J e Eji/T c _ I 

The relations (1III.H IIII.2j) allow one to derive the fictive absorption and stimulated emis- 
sion rates (involved in the "radiative-rate" system) from the known spontaneous emission 
coefficients A,j. 



B. Definition of the bound- free radiative rates 

Assuming that i is a level of an ion with net charge Z* , j with net charge Z* + l, 
photoionization and recombination processes for photon energies between hv and h[y + 5v) 
are described by the kinetic equation 

5 ^t = i^ UuNl ~ T ^ NeN l ~ r j? U » N e N i) (IH.3) 

where r pi , r rr , r srr are the photoionization, radiative recombination and stimulated radiative 
recombination coefficients, which are (supposedly known) functions of the incident wave- 
length. 

When one solves a zero-field collisional-radiative system, one integrates the radiative 
recombination coefficient over the photon energy and get the rate noted as i&. To deal 
with a situation of a Planckian field in equilibrium with the plasma, it would be necessary 
to perform the integration of the photoionization and stimulated recombination coefficients 
in flIII.3j) on the u v Planck distribution, or on the electron Maxwellian distribution, the 
electron energy e being related to the field frequency through + hv = Ej + e. However 
this integration may be avoided since one is interested here not in the determination of real 
photoionization or stimulated recombination rates but in a system involving bound-bound 
and bound-free radiative rates, the solution of which is the Saha-Boltzmann distribution. 
One therefore introduces a unique spectral energy u v computed at threshold v = Eji/h while 
in general one has v > Eji/h. Then, in stationary regime, one changes equation (1IH.3|) into 
the solution at threshold 

(IIL4) 

H 3i 



U v (Eji/h) = —r 



if the Saha-Boltzmann relation (1II.6j) is explicitly used. From the already known recombi- 
nation rate iS, the equations 

fi$ME jt /h) = w 9 -i E . / (III.5) 
R vr 

R-u u (E Tl /h) = &Ej J_ l (III.6) 

provide "fictive" photoionization and stimulated recombination rates respectively. The 
above defined i?? 1 and RJ? are not necessarily accurate values of the real photoionization 
and stimulated recombination coefficients, but they have the two required properties: 

• they depend on the known radiative recombination rate in a simple way 

• the solution of the detailed modified rate equations in this fictive Planckian field 
including bound-bound and bound-free rates as defined in subsections IIII Al and IIIIBI 
is by construction the Saha-Boltzmann solution. 

In the rather usual case where the photoionization cross-sections rapidly decrease at thresh- 
old, the above values are also realistic approximations of the photoionization and stimulated 
recombination rates in a Planckian field. 



IV. A COMPARISON OF VARIOUS PARTIAL RATE EQUATION SYSTEMS 
A. Description of the calculations 

In order to illustrate the proposals of the previous sections, a detailed and configuration- 
average collisional-radiative analysis has been performed in carbon and neon plasmas, in- 
cluding all possible charge states in order to correctly describe a large range of temperatures 
and electron densities. The atomic and collisional computations have been done using the 



HULLAC (Hebrew University Lawrence Livermore Code) suite 27]] . This code includes full 



account of configuration interaction, of particular importance, e.g., in plasmas devised for 



extreme-UV generation 15 



16 



281 ] . It allows for a fast computation of collisional cross- 



sections using the factorization-interpolation method. 

The included configurations in C are the same as in our previous work 
rations in Ne, chosen according to the same criterion, are enumerated in Tabled] This gives 



24|. The configu- 



S 



rise to 150 configurations and 1782 levels in C, where the computation time remain rea- 
sonable while allowing for a valid check of the present criteria for the configuration average 
procedure. In Ne, the 4638-level case give rise to longer calculations which provide a serious 
accuracy check for the proposed criterion. One may estimate that the maximum value of 
the principal quantum number n = 5 may be too low to describe certain quantities, e.g., 
radiative losses, which depend mainly on excited level populations. Several processes such 
as dielectronic recombination may have a large cross-section for high n. However the present 
work is mainly aimed at checking the operation of validity criteria — which of course would 
apply to higher n — and not at providing a reference NLTE computation. Furthermore, a 
test in carbon for N e = 10 16 cm -3 has established that increasing the maximum n from 5 to 
6 changes the average ionization from (Z*J = 3.927 to 3.937 if T e = 10 eV and from 0.934 
to 0.935 if T e — 1 eV, and the computed radiative losses exhibit a less-than-15% change. 
This indicates an acceptable accuracy while keeping the amount of computations — which 
can be considerable when dealing with detailed- level check — moderate. 

The neon computation involves 541 713 radiative deexcitation rates, 1 532 357 collisional 
excitation rates, 2 431 784 collisional ionization rates, 1 144 559 radiative recombination 
rates, and 23 321 autoionization rates. As discussed previously [24|, some transition rates 
are poorly determined in HULLAC. For instance few collisional excitation rates are sensitive 



to inaccuracy in the Sampson fit used 



25 



261 ] . It turns out that, at 30 eV, 6 433 collisional 



excitation rates (0.4%) are negative, 3 402 (0.3%) collisional ionization rates are singular 
24j and 24 (0.1%) autoionization rates are abnormally large. Since these fractions remain 
small, such irregular values may be simply cancelled. 

In order to estimate the influence of such cancellation of some collisional excitation rates 
- for which the higher percentage of unexpected values is observed - - a simple check 
has been performed. It consists in substituting to these unknown rates simple analytical 
expressions relating the collisional deexcitation cross-section to the radiative deexcitation 



rate 



29 



301 ] . For instance, in carbon at 10 eV, 8 190 collisional excitation rates out of 



465 083 (1.8%) are negative when computed with HULLAC. Using a Van Regemorter-type 
formula with an average Gaunt factor equal to 0.2, one may estimate 4 583 additional rates 
(56% of the negative values); the other rates correspond to transitions forbidden by electric- 
dipole selection rules, for which a null cross-section seems to be acceptable. Doing this, the 
average ionization degree at N e = 10 16 cm -3 is (Z*) = 3.9275 while it was 3.9273 when 
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cancelling all negative excitation rates. Therefore, as a rule, the cancellation of few irregular 
rates bears very little consequence. 

A comparison with neon detailed and superconfiguration-averaged values published by 
Hansen et al [l3| is presented in Table HT1 The DLA values also results from a HULLAC-based 
detailed calculation, while MOST is a super configuration code based on a self-consistent field 
(9), [lO]. The "hybrid" model combines the power of statistical average with the accuracy of 



detailed models and is based on the Flexible Atomic Code [311 ] ; contrary to other compu- 
tations, it includes the continuum lowering [14J]. The agreement between these values and 
the present one is satisfactory. As discussed by Hansen et al, the differences arise from the 
diversity in the underlying atomic data rather than from the averaging process. The more 
significant departure with the hybrid model at high electron densities originates certainly in 
the account for continuum lowering. 



B. Using the microreversibility principle in the detailed case 

The solution of the detailed system (III. II) with modified (microreversible) rates as dis- 
cussed in subsection III CI provides a useful accuracy check when comparison is made with 
Saha-Boltzmann. 

This is illustrated for a 10-eV neon plasma in Table HTT1 The average ionization degree is 
computed for the various rate systems detailed at the end of subsection III CI One observes 
that the "collisional" , "radiative", and "complete" rate equations provide a {Z )-value in 
close agreement with Saha-Boltzmann equation, usually agreeing within the nine digits, 
with a small degradation for large electron densities. It is essential to note that this N x N 
(N = 4637) matrix inversion is performed with a remarkable accuracy, the Gauss elimination 
used resulting in a number of floating-point operations scaling as iV 3 . In addition to the (Z* ) 
comparison, the tests on the maximum population difference (e.g., max |iVj(coll) — A^(SB) | 
for 1 < i < 4638) displayed in this table confirms the accuracy of the present rate-equation 
solution. 
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C. Configuration average in the high-temperature regime 



Let us first consider high-T e situations where CA is supposed to be valid. The CA- 
validity test in its different forms has been first applied to a carbon plasma at 10 eV. 
As seen in Table IIV} it turns out that CA is then as a rule acceptable, except for N e > 
10 21 cm -3 where up to 50% divergence is observed. In this case, two of the tests proposed 
here reveal the CA-validity breakdown. While for N e = 10 22 the average dispersion AE a 
is 2 eV, well below the thermal energy, the \Z*\ value computed with "collisional" rates 
(1.053) or "complete" rates (1.070) differ significantly from the Saha-Boltzmann ionization 
(0.718). In this case the ratios (Z* oll ) / (Z$ B ) = 1-47 or / (^sb) = 1-49 are acceptable 

approximations for the CA/detailed ratio (^ca) / (Z*et) = 1-50. However, the "radiative" 
test value (Z r * d ) remain much closer to Saha-Boltzmann, the maximum discrepancy being 
(■^rad) / (^sb) = 0-96 at N c = 10 22 cm -3 . Here the "radiative" test is less sensitive to 



irregularities in autoionization and collisional ionization previously noticed [24| and therefore 
it does not detect accurately the CA-validity breakdown. However a safe way to control the 
CA validity is to compare the charge (Z* n ) (and not (Z* oY ^) to \Z$^) because the latter 
quantity is sensitive to irregularities in all kinds of processes. 

A comparison between detailed and CA results in a 30 eV-neon plasma is presented in 
Table IVl For this temperature the CA approximation is acceptable in the whole density range 
investigated, the error varying from 3% to 1%. However, the energy criterion (1II.4|) . though 
globally satisfied, does not give an accurate picture of the CA validity. For increasing density 
numbers, the CA approximation improves while the criterion (111.41) is less verified. While 
in this table the energy dispersion AE a is evaluated using configuration populations N a 
derived from the CA-collisional-radiative solution, using N a from Saha-Boltzmann equation 
would increase monotonically from 10 -10 eV to 5.7 eV in this density range. 

Nevertheless in the 30-eV neon case, the ratio (Z*n) / (^sb) cannot provide a uniform test 
for CA validity. Noticeably at this temperature, for low N e the most probable ion state at 
thermal equilibrium (very different from the collisional-radiative solution) is the closed-shell 
(Is 2 ) Neix ground state. As a consequence, the ratio (^*n) / (^sb) * s ver y close to unity 
while the CA and detailed charge differ by about 3%. In order to monitor this difference, a 
complementary global information is given by the various central moments of the population 
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distribution 

mt = {{z* - (z*)Y) = J2 m - (z*)Y (iv.i) 

0<k<Z 

where is the ion population of net charge k, and t = 2 or 3 here. The corresponding 
results for a 30-eV neon plasma are presented in Tables [VTI and IVlTI One may check that, 
at N e = 10 12 cm -3 for instance the ratio m2(coll)/m2(SB) is 0.93 (the ratio formed with 
"radiative" or "complete" rri2 being very similar) while the m2(CA)/m2(det) is 0.86. Even 
(much) stronger discrepancies could be obtained using third central moments. Therefore the 
use of second (and possibly third) central moments together with the average charge and 
comparison of "collisional" or "complete" data to Saha-Boltzmann data provide as a rule 
an efficient and uniform test for the CA validity. 



D. Configuration average in the low-temperature case 

In figure [1] are presented the average charge (Z*} and the second central moment ni2 
( 1IV.1I) for a 10-eV neon plasma, both in the detailed-level scheme and in the CA approxima- 
tion. The figure includes Saha-Boltzmann results too. For densities large enough (at least 
10 20 electrons/cm 3 ), the detailed data obtained with the collisional-radiative code (crosses) 
are identical to the Saha-Boltzmann values (circles), indicating that LTE is then reached. 
One notes that Saha-Boltzmann data are almost identical in the detailed (circles) and the 
CA case (dashed line). This indicates that the averaging on energies then performed in 
the Saha-Boltzmann equation bears very little influence on populations. Conversely, the 
CA (crosses) and detailed (solid line) results obtained with the collisional-radiative system 
are significantly different, indicating that the averaging procedure on rate equations is not 
of little consequence. As expected the CA validity worsens at low T e , the difference being 
(Zq^ — (Z* et ) = 0.92 at A c = 10 20 cm" 3 . As seen in the lower part of this figure, the 
second moment m 2 exhibits similar properties, e.g., concerning the LTE convergence in the 
detailed case and the collisional-radiative vs Saha-Boltzmann differences. One can easily 
check that, if few charge states are populated, when \Z*) is integer (resp. half-integer), 
m>2 is close to a minimum (resp. maximum). This explains the oscillations observed on this 
figure. 

In order to check the breakdown of CA validity without performing any detailed-level 
computation, a first criterion is provided by the average-energy dispersion flll.4j) . The first 
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member of this equation is plotted as a function of N e in figure [2j It appears that the energy 
dispersion (AE) is here rather large, in the 2-5 eV range, to be compared to T e = 10 eV. 
Furthermore, the dispersion calculated with the collisional-radiative system solution is al- 
ways greater than 1.9 eV, and peaks at 5.2 eV for N c = 10 20 cm -3 which is precisely the 
maximum of the CA versus detailed difference on \Z*Y If one estimates that the criterion 
(1II.4j) is not fulfilled for (AE) = 2.6 eV, it has the expected behavior. The energy dispersion 
(AE) may also be computed with populations from the Saha-Boltzmann equation or from 
the modified rate equations proposed in subsection III CI As seen in figure [2, it is smaller than 
the dispersion from collisional-radiative solution (except in the N e ~ 10 15 -10 16 cm -3 range) 
and particularly for electron densities well below 10 13 cm~ 3 where this quantity tends to 
zero. This can be explained as follows, e.g., for the "collisional" rate equations. Ignoring the 
collisional excitation and de-excitation processes that do not change the ionization degree, 
the ionization balance depends on the collisional ionization rate from the charge state z to 
2 + 1, R( C1 > N z N e , and on the three-body recombination process from z+1 to z, R^N z+1 Nl 
At equilibrium, this gives a population ratio N z /N z+ % ~ N e R( 3br ^ / R( a \ vanishing at low N e . 
Therefore the highest charges are then favored: it should be 10 at very low N c , but since 
the closed-shell Is 2 Ne IX is very stable, one gets \Z*) = 8 over a large range of densities. 
A similar analysis may be performed for the "radiative" equilibrium where photoionization 
and radiative recombination are included. In these cases, since the mainly populated con- 
figuration is Is 2 with AE a = 0, the resulting average energy dispersion is very small. In 
the contrary, in the collisional-radiative case, the low-density balance (coronal plasmas) is 
governed by collisional ionization and radiative recombination and then one gets a constant 
nonzero N z /N z+ %, and therefore (Z*) ^ 8. For this rate system, the average dispersion 
(AE) is nonzero, as seen in figure [2j 

To check how the proposed criterion based on the solution of a partial rate equations be- 
haves, one may consider the variation of (Z* h ) — (Zg B ) displayed in figure [3j It is expected 
that, when this difference is large, the CA approximation will fail. It first appears that the 
"radiative" solution as proposed in section II III remain rather close to the Saha-Boltzmann 
solution, and cannot in this case provide a test for the CA validity. However, the "complete" 
solution involving all the microreversible rates (collisional, autoionization and radiative in a 
fictive Planck field) presents a maximum difference on (Z* n ) — (Z§ B ) at iV c = 3 x 10 20 cm -3 , 
precisely where the CA-detailed departure is maximum. Furthermore, this (Z* n ) — (Z§ B ) 
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value, close to 1, provides then a reasonable quantitative estimate of the CA-detailed dif- 
ference, at variance with the simple criterion (111.4j) . The drawback of the present criterion 
lies here once again the the low-density region, where all the computed {Z*} are close to 8, 
and therefore the differences plotted on figure [3] tend to zero. In this particular case, one 
must resort to the usual energy criterion (III.4p . This is an opposite situation to the carbon 
case j^J]: for instance, if N e = 10 12 cm -3 and T c = 1 eV, while the energy criterion (111.41) is 
fulfilled by large, the average charge difference \Z* oll } — (Z| B ) is close to 1, indicating that 
CA approximation does fail. 

An alternate analysis may be performed on the second or third central moments computed 
with the solution with rates as defined in section HTT1 and compared to the analogous moments 
from Saha-Boltzmann solution. It has been checked that the conclusions remain similar, with 
a maximum difference around N e = 10 21 cm -3 . Again the "radiative" solution remain closer 
to Saha-Boltzmann than the "collisional" or "full" solution. 



V. CONCLUSION 



It has been demonstrated that the configuration-average validity for collisional-radiative 
rates, which leads to considerable simplification of these equation systems, may be controlled 
by comparison of several variants of a modified rate system with Saha-Boltzmann solution. 
These tests have been performed on the average ionization degree as well as the second 
and third central moments. The radiative test appears to be easier to verify than the colli- 
sional or complete tests: This demonstrates that the CA breakdown is induced mainly by 
large dispersions in collisional rates inside a given pair of configurations and not in radiative 
rates. Tests have been performed in carbon and neon, where atomic and collisional data 
were provided by the HULLAC code. Going from C to Ne amounts to increase the number 
of levels from 1782 to 4638 (and from 150 to 291 configurations) without impairing the effi- 
ciency of the proposed test. A detailed analysis in a carbon plasma has proven that failure 
of the CA at T e = 1 eV may be correctly detected by the test based on the modified-rate 
system, while failure of CA in a 10 eV-neon plasma at low density is better diagnosed by 
the intuitive criterion based on the average energy dispersion inside a configuration. This 
situation prevails when the modified-rate equations solution is dominated by closed-shell 
ions such as Neix. Therefore, both criteria appear quite complementary. Nevertheless, the 
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criterion based on the average charge difference may provide a semi-quantitative estimate of 
the CA-versus-detailed charge difference, at variance with the dispersion energy criterion. 
Further developments include the use of such systems in the determination of physically im- 
portant quantities such as radiative losses, the opacity, and the emission of non-equilibrium 



plasmas 



32 



331 ] . It is interesting to check how the various forms of the present criterion 



can control the relevance of CA approximation when dealing with such radiative properties. 
As a possible example of application, a configuration averaged code based on the HULLAC 
suite has been recently applied in our group to the EU V emission of xenon plasmas [34 ] . 
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APPENDIX A: VALIDITY OF THE CONFIGURATION AVERAGE: AN EX- 
AMPLE 

In the special case where the rates inside every pair of configurations are proportional to 
the final level degeneracy 

R ij = R a p— ViEa,j£p, (A.l) 

9/3 

the solution of the detailed-level rate equations can be straightforwardly derived from the 
solution of the configuration-average rate equations. Let us assume that the populations N a 
are the solutions of the CA rate equation fill. 31) . If we define level populations according to 

Ni = N a —, Nj = N/3— tea }3 e(3 (A.2) 

9a 9/3 

then it is easy to show that such populations obey the detailed rate equation (III. II) . One 
has according to (lA.lj) and (1A.2I) 



V {-RijNi + RjiNj) = Y] (-R aP ^N a + R Pa ^N^\ (A.3) 



J 3 
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where a (resp. (3) is the configuration containing i (resp. j). In the above sum, when the 
level j belongs to the same configuration as i (fl — a), the j-term vanishes identically; when 
j belongs to another configuration (/3 ^ a), the sum over j ^ a may be written J2/3^ a 
and then, since YljepOj/Qp = 1' one § e ^ s after performing the j-sum 

{-RiA + RjiNj) =-J2 (~ R ^ a + RpvNp) = (A.4) 

3 9a 

because of the A^ a , Np satisfy the CA rate equations (111.31) . In a case where ( 1A.1I) is fulfilled 
the microreversibility-based (or "thermodynamic") test checking the CA validity proposed 
in this paper is satisfied: the populations obtained with partial rate equations are identical 
to Saha-Boltzmann. 



APPENDIX B: AVERAGE OF MICRO REVERSIBLE RATES WHEN THE EN- 
ERGY DISPERSION IS SMALL VERSUS THE ELECTRON THERMAL ENERGY 

In this appendix it is proved that, if the energy dispersion inside every configuration is 
much smaller than the electron thermal energy ( 111.71) , then the microreversibility condition on 
detailed levels implies microreversibility on configuration- averaged rates too. Let us assume 
that the processes p and q (e.g., collisional ionization and three-body recombination) obey 
the microreversibility condition 

R^Nf B =R^Nf B (B.l) 

for every i e a, j G /3, where the populations Nf B , Nj B obey the Saha-Boltzmann equation 
( 1II.6|) . Then, using the average rate definition (1II.2j) one gets 



_ m (B.2) 

= ^£fte ( - Ej ' /T - > < ) 

and if it is possible to identify the energy difference Eji with the difference of the average 
configuration energy 

\Eji — Ep a \ <§C T e (B.3) 
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then after substitution of Ep a to Eji, the rate expression (1B.2j) involves the Saha-Boltzmann 
ratio of the configuration populations 

ie/3 (B.4) 

atSB 

- di^s-E^/nM _ l 2l p(g) 

where the Saha-Boltzmann populations iV^ B computed with average energies have been 
introduced. This means that, in this case, the microreversibility condition holds for configu- 
rations too. However, conditions fill. 71) or (IB.3I) assume that the maximum energy dispersion 
on ewer?/ pair of configurations must be smaller than T e , which is a much stronger condition 
than ()II.4j) which only requires the average energy dispersion to be much less than T e . 
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TABLES 



TABLE I: Number of configurations and levels considered for the collisional-radiative calculations in 
the neon plasma. For Nei to Nevin, the included configurations are ls 2 {2s2p} k , ls 2 {2s2p} k ~ 1 Nl, 
where {2s2p}i stands for 2s 2 2p ? ~ 2 , 2s2p ? ~ 1 , 2p ? ~ 2 . The configurations considered in Neix are 
Is 2 and IsNl, those in Nex are Nl. One has N < 5, / < N - 1 for each ion. The Ci-Cvii 
computations include configurations isoelectronic to Nev-Nexi respectively. 



Ion 


Nei 


Neil 


Nein 


Neiv 


Nev 


Nevi 


NeVll 


Nevin 


Neix 


Nex 


Nexi 


Total 


Configurations 


25 


38 


39 


39 


39 


39 


27 


14 


15 


15 


1 


291 


Levels 


157 


501 


994 


1204 


1004 


513 


166 


24 


49 


25 


1 


4638 



TABLE II: Average ionization level in a neon plasma at T e = 25 eV and 50 eV. DLA is a detailed 



model, while MOST is based upon super configuration averaging 13 1. The code "hybrid" is a 



partially detailed and partially averaged model based on FAC code 
are the present results, in detailed and configuration average form. 
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311 ] . The last two columns 



T 


iV c (cm- 3 ) 


DLA 


MOST 


hybrid 


detail 


CA 


25 


10 16 


5.247 


5.87 


5.227 


5.530 


5.712 




10 18 


6.050 


6.17 


6.101 


6.352 


6.406 




10 20 


6.158 


6.23 


6.291 


6.183 


6.303 




10 22 


3.403 


3.78 


4.200 


3.083 


3.248 


50 


10 16 


7.625 


7.45 


7.618 


7.604 


7.607 




10 18 


7.739 


7.80 


7.797 


7.791 


7.791 




10 20 


7.901 


7.90 


7.917 


7.905 


7.905 




10 22 


6.275 


6.31 


6.706 


6.134 


6.194 
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TABLE III: Accuracy check of the collisional-radiative (CR) detailed-level solution for neon at T c = 10 eV. The average charge (Z*) is 
calculated using three kinds of microreversible rates: collisional rates ((Z* on )), radiative rates including absorption and stimulated emission 
in a Planckian field ((Z* ad )), all these rates plus autoionization and dielectronic recombination rates ({Z* n )). (Z$ B ) is the average charge 
at thermal equilibrium as derived from Saha-Boltzmann and should be equal to the previous (Z*) values at infinite numerical accuracy. 
The additional test <5 ra( j (resp. 5 co \, <5 a n,) is the maximum difference on the ion-level populations between the "radiative", (resp "collisional", 
"complete") system and the Saha-Boltzmann solution. 



JV e (cm -3 ) 




<^oii> 


<*a*»> 


<^S*B> 


<5rad 




<5 C ol 




Sail 




10 12 


7.712640377 


7.712640378 


7.712640377 


7.712640377 


2.0 x 10~ 


-10 


2.7 x 10- 


-9 


1.5 x 10" 


■10 


10 14 


6.639625837 


6.639625830 


6.639625838 


6.639625838 


1.9 x 10~ 


-10 


3.4 x 10 


-9 


4.6 x 10- 


•li 


10 16 


5.491015460 


5.491015458 


5.491015465 


5.491015460 


1.9 x 10- 


-12 


7.2 x 10~ 


•10 


1.3 x 10" 


-9 


10 18 


4.260724652 


4.260724657 


4.260724656 


4.260724652 


6.5 x 10- 


-13 


1.5 x 10" 


-9 


1.3 x 10" 


-9 


10 20 


2.847355968 


2.847355980 


2.847355979 


2.847355968 


2.1 x 10' 


-14 


1.1 x 10" 


-9 


1.2 x 10" 


-9 


10 22 


0.719532736 


0.719535868 


0.719535378 


0.719532736 


4.1 x 10- 


-14 


6.9 x 10- 


-7 


7.0 x 10" 


-7 



TABLE IV: Average ionization stage in a 10 eV-carbon plasma as a function of the electron density 
N c . The last two columns are the configuration- average (CA) charge and the detailed- level charge 
derived from the collisional-radiative (CR) equations. The discrepancy between these values is 
correlated to the discrepancy between the Saha-Boltzmann charge (column 5) and the average 
charge derived from various modified rate-equation systems: collisional rates (column 2), radiative 
rates including a Active Planckian field (column 3), all these rates plus autoionization (column 4). 
Data in columns 2-5 are obtained in the CA framework. 







Modified-rate system 




CR system 


(cm -3 ) 


<^oli> 


( Z rad) 


<35i> 


<^S*B> 


(Z*ca) 


(zL) 


10 12 


4.000003693 


4.000003671 


4.000003671 


4.000004022 


3.7363 


3.7441 


10 14 


3.999998275 


3.999998275 


3.999998275 


3.999998278 


3.7563 


3.7616 


10 16 


3.999823807 


3.999823806 


3.999823806 


3.999823805 


3.9273 


3.9273 


10 18 


3.982626435 


3.982615120 


3.982627955 


3.982613441 


3.9679 


3.9679 


1Q 20 


3.225872703 


3.197765566 


3.230477092 


3.194905628 


3.2282 


3.1902 


10 22 


1.052601006 


0.689418858 


1.069768540 


0.717721378 


1.0698 


0.7113 
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TABLE V: Average ionization degree in a neon plasma at T e = 30 eV with respect to the electron 
density A c . Column 2 involves detailed levels, while columns 3-7 contain charges which are com- 
puted in the configuration average scheme. Columns 2 and 3 are the collisional-radiative solution, 
while columns 4-6 are variants of "microreversibility check", based on a collisional, radiative, or 
complete set of rates respectively. The closer these values are with respect to the Saha-Boltzmann 
(■£*) value (column 7), the better is assumed to be the configuration average. Column 8 is the aver- 
age energy dispersion of levels inside a configuration with a ponderation by the collisional-radiative 
configuration populations. 



A c 


CR syst 


em 


Modified rate sysl 


;em 


Saha- 


(A£ CR ) 


(cm -3 ) 


Detailed 


CA 


Collisional 


Radiative 


Complete 


Boltzmann 


(eV) 


10 12 


5.7919 


5.6236 


8.000009396 


8.000009460 


8.000009460 


8.000010081 


0.069 


10 14 


5.7911 


5.6540 


7.999998415 


7.999998416 


7.999998416 


7.999998422 


0.164 


10 16 


6.2120 


6.3134 


7.999832154 


7.999832154 


7.999832154 


7.999832154 


2.18 


10 18 


6.8256 


6.8456 


7.983410902 


7.983409743 


7.983409658 


7.983407736 


1.13 


10 20 


6.9131 


6.9400 


7.180274366 


7.172115334 


7.180822653 


7.168028134 


0.97 


10 22 


3.6140 


3.6603 


4.471505865 


4.090408345 


4.436454002 


4.047538516 


5.92 



TABLE VI: Second central moment ({Z* — (Z*)) 2 ) in a neon plasma at T e = 30 eV as a function 
of N e . Refer to Table M for details. The average energy dispersion {AEqr) is not repeated. 



A c 


CR system 


Modified rate system 


Saha- 


cm -3 ) 


Detailed 


CA 


Collisional 


Radiative 


Complete 


Boltzmann 


10 12 


0.6128 


0.5283 


0.000009430 


0.000009494 


0.000009494 


0.000010114 


10 14 


0.6137 


0.5433 


0.000001773 


0.000001773 


0.000001773 


0.000001780 


10 16 


0.5043 


0.4528 


0.000167828 


0.000167828 


0.000167828 


0.000167829 


10 18 


0.3288 


0.3137 


0.016394148 


0.016396461 


0.016396557 


0.016400426 


10 20 


0.4700 


0.4377 


0.442653567 


0.456411086 


0.440578198 


0.459907235 


10 22 


0.7707 


0.7484 


0.785410789 


0.787926179 


0.811753013 


0.823178341 
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TABLE VII: Third central moment {{Z* - (Z*)) 3 ) in a neon plasma at T e = 30 eV as a function 



of N e . Refer to Table E] for details. 



N c 


CR system 


Modified rate syst< 


am 


Saha- 


( CTT1 ~~ ^1 


Detailed 


CA 


Collisional 


Radiative 


Comnlete 


Boltzmann 


10 12 


0.0647 


0.2229 


0.000009396 


0.000009460 


000009460 


000010080 


10 14 


0.0685 


0.2146 


-0.000001585 


-0.000001584 


-0.000001584 


-0.000001578 


10 16 


-0.0500 


-0.0416 


-0.000167786 


-0.000167786 


-0.000167787 


-0.000167787 


10 18 


-0.0289 


-0.0212 


-0.016009480 


-0.016014093 


-0.016014161 


-0.016021878 


10 20 


-0.0844 


-0.0508 


-0.097745732 


-0.118986694 


-0.093018605 


-0.119420248 


10 22 


0.2376 


0.3694 


-0.125419985 


-0.024429621 


-0.109647319 


-0.010892532 
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FIGURES 



FIG. 1: Average ionization (Z*) and second central moment 777,2 = — {Z*)) 2 ) for a 10-eV 
neon plasma, in the detailed level scheme (symbols) and in configuration average (lines). Popula- 
tions are computed either with the collisional-radiative code or with the Saha-Boltzmann equation. 
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FIG. 2: Average configuration energy dispersion (AE) = ^ Q N a AE a in a T c = 10 eV neon plasma. 
AE a is the energy dispersion inside the configuration a, and N a is the a-population calculated with 
various rate-equation systems and with Saha-Boltzmann equation. The smaller these dispersions 
are versus T e , the better is assumed to be the configuration average approximation. 
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FIG. 3: Difference between the average charge calculated with various rate-equation systems and 
the average charge obtained from the Saha-Boltzmann equation. Calculations are performed in a 
10-eV neon plasma in the configuration average scheme. 
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